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I ■ We investigate the recent CDF measurement of the Bt(Bs —^ 1^^ 1^) which shows excess over 

r^. the Standard Model. We consider minimal supergravity motivated models (mSUGRA)/CMSSM 

^ ! and grand unified models, SU(5) and SO(IO). In the grand unified models, the neutrino mixings 

provide an additional source of squark fiavor violation through the quark-lepton unification. In 
the context of minimal SU(5) model, we find that the new CDF measurement has imposed a 
Q ■ lower bound on the branching ratio of r — )■ yU7 for a large CP phase in the Bg-Bg mixing. Recall 

CN . that there have been indication for a large CP phase in Bg mixing from Bg — )■ J/ip4' (Tevatron 

r~ ! and LHCb) and dimuon asymmetry (DO). We also predict Br(r — t- fir]) for the possible range 



of values of Br(r — )■ 7x7). 



1 Introduction 

Recently, the CDF collaboration has reported the measurement of the rare decay process of Bg 
meson to fi~^fi~ [T] , and the branching ratio of the decay is 

Bt{B, ^ /^+/i~)'^°^ = (1-8+^]) X 10-^ (1) 

The 90% confidence level (CL) range is 

4.6 X 10"^ < Br(5, -^ /iV")^°^ < 3.9 x 10"^ (2) 

The measured branching ratio at the 90% CL is deviated from the standard model (SM) 
prediction [21 13] : 

Bi{Bs -^ f^^f-i'f^ = (3.19 ± 0.19) X 10"^. (3) 

The SM prediction above is calculated from the ratio of the branching ratios and the Bg mass 
difference, AMg. The discrepancy cannot be explained by the hadronic uncertainties, which 
comes from the uncertainty on the bag parameter Bb^ = 1.33 ±0.06 [3j. Thus, this discrepancy 
implies the existence of a new physics (NP) beyond the standard model. The excess can be 
soon verified at the LHC |1]. 

The excess of the Bs — )• fi'^fi~ rare decay can be reproduced in many NP models. Su- 
persymmetry (SUSY) is one of the promising candidates for the new physics. In the minimal 
SUSY extension of the standard model (MSSM), the rare decay of Bs — )■ fi^fi~ is induced 
by the neutral Higgs mediated flavor changing operator [5]. The operator is proportional to 
ta.n^ f3/m\, where tan/3 is the ratio of vacuum expectation values of the up-type and down- 
type Higgs fields, and rriA is the mass of the CP odd Higgs field. Therefore, the branching 
ratio is proportional to tan^ /3, and the excess of this rare decay process implies a large value of 
tan P for the experimentally allowed parameter region of the model, where the masses of SUSY 
particles (especially for stops) are bounded from below. The flavor changing Higgs coupling is 
induced by finite corrections in the down-type quark mass matrix, and it can be generated even 
in the minimal flavor violating models where the CKM (Cabibbo-Kobayashi-Maskawa) quark 
mixing is the only source of the flavor violation. In this case, the operator is generated by a 
chargino-stop loop. 

In the SUSY grand unified theories (GUTs), flavor mixing in the lepton sector (i.e. neutrino 
mixing) can cause a squark flavor violation due to the quark and lepton unification [6| 17]. The 
squark flavor violation can induce an additional contribution to the Higgs mediated flavor 
changing coupling via a gluino-sbottom loop. In the minimal type of SU(5) GUT, the right- 
handed down- type quark and the left-handed lepton doublet are unified in a multiplet, the 



right-handed down-type squarks have the flavor violating source [8], and the left-handed squarks 
(i.e. both up- and down-type) can have the flavor violating sources in SO(IO) GUTs, where 
the size of the squark/slepton flavor violation can be related to the enhancement of the proton 
life time [9] . The detail investigation of the flavor violating processes are important to find the 
footprint of the SUSY GUTs [TOl [IH [121 [I3] . 

The right-handed squark flavor violation with a large tan /3 can induce a sizable Bs-Bg mixing 
amplitude by a flavor changing Higgs mediation rather than the box diagram contribution 
P^ [TS| [TB| [T7] . The h-s flavor violation, which is suggested by B decay experiments (the CP 
violation in B^ -^ J /'4"P decay [HI [19] and dimuon asymmetry in semileptonic B decays [20] ) , 
can be related to the r-/i flavor violation, such as r — )■ /i7, r ^i- 3fi and r — )■ /ir/, in SUSY 
GUTs. Due to the bounds of the flavor violating r decays, especially for r -^ fi'j, the size of 
the b-s flavor violation and the CP odd Higgs mass can be bounded to obtain a large CP phase. 
As a result, the existence of a large CP phase for Bg-Bg mixing can constrain Bt{Bs — ?■ fJ'^fJ'^) 
in the case where the quark and lepton unification is manifested. In fact, we obtain that 
the Bt{Bs — ;■ fi^fJ^') should be larger than O(10~®) [21J. Inversely speaking, the branching 
ratio measured by CDF gives a bound on the flavor violating r decays. In the paper, we will 
investigate the bound of r decays in SUSY GUTs. 

This paper is organized as follows: in section II, we discuss Bg — t- fJ^^fJ^" in minimal super- 
gravity (mSUGRA), in section III, we discuss large Bg mixing phase in SUSY GUT models, in 
section IV, we discuss the flavor violating r decays and we conclude in section V. 

2 Bs ^ fi^fi~ in minimal supergravity (mSUGRA) 

The rare decay of Bg — )■ ^'^ ^~ can be generated even if the CKM quark mixing is the only source 
of flavor violation in MSSM. We first study the parameter space in minimal supergravity model 
(mSUGRA) (constrained MSSM (CMSSM)) [22] to realize the CDF measured Bg -^ /i+^-. 

In mSUGRA, the parameters are the universal scalar mass (mo), the unified gaugino mass 
{fTti/2)i the universal trilinear scalar coupling (Aq), the ratio of the Higgs vev (tan/3), and the 
sign of the Higgs mixing parameter yU. The current experimental constraints, 6 — )■ 57 and the 
lightest Higgs mass bounds raise the squark masses. As a result, lower tan/3 is disfavored to 
obtain Br(5, -> ^1+ ^i') > 10"*^ [23j. 

Recent results from ATLAS with 165 pb^^ of luminosity show that the maximum ruled out 
value of mi/2 is ~ 450 GeV [2^ , which corresponds to squark, gluino masses ~ 1 TeV. 

In FigUl we show a typical parameter space for mSUGRA, with tan/3 =40 and 50. The 
model parameters are already significantly constrained by different experimental results. The 
most important constraints for limiting the parameter space are: (i) the light Higgs mass bound 
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Figure 1: The mSUGRA/CMSSM parameter space is shown. The CDF allowed range of 
BR(i?s — )• fi^fi^) is shown by the green shaded region. The other shaded regions and lines are 
described in the text. 



of rufio > 114.4 GeV from LEP [23] (red dotted line shows the contours of ruho = 114.4 GeV, 
as calculated using FeynHiggs-2 . 6 . 5 [2S]); (ii) the b ^ sj branching ratio [27J (95% CL 
excluded in the yellow shaded region in Fig. [1]); (iii) the 2a bound on the dark matter relic 
density: 0.106 < ^cBMh'^ < 0.121 from WMAP [28j (blue region in Fig. H]); (iv) the bound 
on the lightest chargino mass of m-± > 103.5 GeV from LEP [29| (region left to the black 
dashed line is excluded) and (v) the muon magnetic moment anomaly a^ = {g^ — 2)/2 (pink 
shaded region in Fig. [1] is within 2a of where one gets a 4(j deviation from the SM as suggested 
by the experimental results [30] and the analysis in [3T]). We also show the 2a contours 
from 3.2a deviation based on [32] by slanted dashed purple lines. These two references use 
recent changes in the hadronic contribution to calculate the leading order hadronic contribution. 
Assuming that the future data confirms the a^ anomaly, the combined effects of (7^ — 2 and 
m-± > 103.5 GeV then only allows ^ > 0. The grey shaded region in Fig. [T] is excluded 
for not satisfying the electroweak symmetry breaking condition, while the brick-colored region 
is excluded because the stau is lighter than the neutralino hence neutralino cannot be the 
dark matter candidate. The red solid line shows the neutralino-proton elastic scattering cross- 
section contour of 9 x 10^^ pb which approximately is the bound from XENONIOO [33] for 
70 GeV neutralino mass. Note, however, that the theoretical cross-section can easily have 
large uncertainties due to the hadronic factor determinations, the dark matter profile and the 
galactic velocity distribution. We show the CDF 90% CL contour of Bt{Bs — )■ fi^fi'), as in 
Eq.(|2]), as a green shaded region. The maximum allowed values of mi/2 and tuq go up to 



Table 1: Sparticle mass ranges for tan/3=50 and 40 for 90%CL and la allowed values of the 
Bt{Bs — 7- fJ^^fJ^") in the dark matter allowed regions. The lower bound for tan/3 = 40 is fixed 
by the 6 — )■ S7 constraint, hence the same for both 90% CL and Icr. 
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0.46 
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2917 
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0.52 
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1952 
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1.05 
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1593 


1804 


595 


464 
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891 


(la) 


0.58 


1341 


1143 


1141 


913 


1042 


351 


252 


246 


474 



~ 1100 GeV for Aq = and tan/3 = 40. This covers the whole neutralino-stau coannihilation 
band favored by dark matter. If we increase tan/3 to 50, the region allowed by the 90% CL 
range of Bt{Bs — )■ yU^/U^) increases and the maximum allowed value of mi/2, that satisfies the 
dark matter constraint, goes up to 1.6 TeV. In addition, we also show the la contour, Eq.([T]), 
as dashed light-green lines in both figures. The lowest tan/3 value, for Aq = 0, allowed by the 
limits on the Higgs mass and the b ^ s'j branching ratio, and within the CDF la region of 
Bt{Bs — 7- yU^/i^) is about 30, but using the 90% CL range tan/3 can be as low as about 20. 

In Table [1] we list the upper and lower spectrums of the sparticle masses in the mSUGRA 
model for tan/3=50 and 40 for the 90%CL and la allowed values of the Bt[Bs — > fi^fi^), when 
we satisfy the dark matter and other constraints. 

Because Bt{Bs — > Ai+/i~) depends on the CP odd Higgs mass niA and the Higgsino mass fi, 
the constraints of the SUSY parameters from the CDF measurement of the branching ratio is 
different in the case of the non-universal Higgs mass (NUHM) boundary condition. The recent 
analysis of Bt{Bs — )■ fi^fi^) in NUHM can be found in [31] (See also [35] for earlier analysis.). 



3 Large B, mixing in SUSY GUTs 

In SUSY GUT theories, it is often assumed that the SUSY breaking sfermion masses are flavor- 
universal, but the off-diagonal elements of the mass matrices are generated by the loop effects. 
The FCNC sources are the Dirac/Majorana neutrino Yukawa couplings, which are responsible 
for the large neutrino mixings ^. Since the left-handed leptons (L) and the right-handed 
down- type quarks (D^) are unifled in 5, the Dirac neutrino Yukawa couplings can be written as 
Y^yijbiNjH^, where A^^ is the right-handed neutrino. The flavor non-universality of the SUSY 



breaking D'^ masses is generated by the colored Higgs and the N'^ loop diagram [8], and the 
non-universal part of the mass matrix is ^M"^^ c^ — g^j (S^q + Al)Yi,Yj} \n{M^/MHc)y where M^: 
is a cut-off scale (e.g. the Planck scale), Mh(^ is a colored Higgs mass, rriQ is the universal scalar 
mass and Aq is the universal scalar trilinear coupling. The left-handed Majorana neutrino 
coupling LLAl (A^ is an SU(2)l triplet) can also provide contributions to the light neutrino 
mass (type II seesaw [36]), and can generate the FCNC in the sfermion masses when the 
fermions are unified. 

As a convention in this paper, we will call the model with the FCNC source arising from the 
Dirac neutrino Yukawa coupling as the minimal type of SU(5). In this case, the off-diagonal 
elements of 10 {Q, If^, E^) representations are small because they originate from the CKM 
mixings. In a competitive model which we call the minimal type of SO(IO), the Majorana 
couplings, which contribute to the neutrino mass, generate the off-diagonal elements for all 



sfermion species since the Majorana couplings fijLiLjA^ can be unified to the fijlGi 16 j 126 
coupling [37]. The detail can be found in Ref . [TT| [2T| . 

The NP contribution of the Bg-Bs mixing amplitude can be described as 

M!2 = Mf2,sM + Mi^2,NP = Cs Mi^2,sM e'^^^^ , (4) 

where Cg is a real positive number. From the measurement of the mass difference, AMg = 
2|M^2l) the experimental result is consistent with C^ = 1. Even if C^ = 1, there is room for 
new physics because of the phase freedom, (pB^ ■ The phase can be measured by Bg — )■ J/ipcj) 
decay and the dimuon asymmetry from semileptonic B decays. In the SM, the CP violation of 
Bg ^ J/ixp is tiny (2/3f^^ ~ 0.04). The CP phase 2/3,(= 2(/3f^ + 0bJ) and the decay width 
difference AF^ have been measured at the Tevatron [18] and the LHCb [19], and a large phase 
is allowed. The dimuon asymmetry reported by DO implies a large CP phase in Bg-Bg mixing 
[20]. 

By definition in the Eq.Q, we obtain 

sm 0B, = ^fj ' (5) 

where A^^ = |M^*2,npI and A^^ = |M;f2,SMl- ^^ ^^e case of Cg ~ 1, we obtain 2sin0B^ ~ 
A^^ /A^^. The large dimuon asymmetry reported by DO requires A^^ /Af^ ~ 0(1). 

In SUSY GUTs, where quark and lepton unification is manifested, such large values of 
A^^ /A^^ indicate a large lepton flavor violation such as r — )■ /i7. Therefore, in order to satisfy 
the current experimental bound 



Br(r -^ fx-f) < 4.4 x 10-^ (6) 
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Figure 2: We plot the correlation of Br(r — )■ jj^y) and Bi{Bs — )■ /i+yU ) for A^^ /Af^ 
the case of the minimal type of SU(5) for universal boundary condition. 
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the parameters in the model are constrained. Since Br(r — )■ fij) is proportional to tan^/3, the 
constraint is more severe for a larger tan/3. The CDF reported Bt{Bs — )■ fi^fj^^) actually prefers 
a large tan/3. 

In the MSSM, the Bg-Bg mixing amplitude is induced by the box diagram contribution and 
the Higgs mediated contribution. The box contribution does not depend on tan/3 explicitly, 
while the Higgs mediated contribution is proportional to tan^ (5/m\. Therefore, for a large 
tan/?, the Higgs mediated contribution can dominate over the box contribution, and it can 
induce a sizable value of A^'^ /A^^ satisfying the experimental bound of Br(r — )■ /i7) in SUSY 
GUTs. The flavor violating Higgs coupling can induce not only a large contribution to the 
Bg-Bg mixing amplitude but also to Bg — )■ /U''"/i~. As a result, for a given size of A^^ /A^^, 
Br(r — !■ /i7) and Bi^Bg — )• fi^fi~) are bounded. 

In FiglJl we show the correlation of Br(r — )■ fi'-f) and Br(i?s — )■ fj^^fi^) for A^^ /A^^ = 0.5 in 
the case of the minimal type of SU(5). In the plot, the Higgsino mass fi and the CP odd Higgs 
mass rriA are varied assuming that the SUSY breaking Higgs masses mn^ and rriH^ are different 
from the universal sfermion mass (mo = 1 TeV). We choose the gaugino mass mi/2 = 500 GeV 
and tan/3 = 40. For a given value of Br(r — )■ fi'y), BT{Bg — )■ fi~^fi~) is bounded from below, and 
vice versa. In this example, in order to satisfy the CDF result Bi^Bg — > fi^ ^~) < 3.5 x 10~^, 
the branching ratio of r — )■ /i7 has to be larger than 10~^. Bi{Bg — )■ ^^ ^~) is also bounded 
from below due to the Br(r — )■ /i7) constraint. 
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Figure 3: The boundaries of Br(r — )■ fij) and Bt{Bs — )■ /i+/i ) correlation regions are shown 
for different input parameters uiq and mi/2- We choose A^^ /A^^ = 0.5 (left plot) and 1 (right 
plot). 

In Figj3l we plot the boundaries of Br(r — )■ fi'-f) and Bt{Bs — )■ fi^fJ'^) correlation regions in 
the same manner as the plot in Fig 121 while changing the input parameters mo and mi/2- We 
plot two cases A^^ /Af^ = 0.5 and 1. Surely, a larger A^^ /A^^ = 1 provides more stringent 
constraints. The squark mass is less dependent on mo due to the gluino loop compared to the 
slepton mass. Therefore, larger rrio gives a smaller Br(r -^ fi'-f). In the case of A^^ /Af^ = 1, 
mo = 500 GeV is already excluded by By{Bs — ?■ fi^fJ'^) and Br(r — )• fi'y) bounds. 

These bounds are significant when the quark and lepton unification is manifested. In the 
minimal type of SU(5) GUT where the squark flavor violation is generated from Dirac neutrino 
Yukawa coupling, the unification is manifested. In the case of SO(IO) GUT, the bounds can 
be relaxed by a choice of symmetry breaking vacua [9]. Therefore, the accurate measurements 
of By{Bs — )■ fi^fi") and the Bg-Bg mixing phase are informative to distinguish the symmetry 
breaking vacua. 



4 Flavor violating r decays 

The rare decay process Bg — )■ fj^'^fi~ is dominantly generated by the Higgs mediated diagram. 
The leptonic version of the Higgs mediated diagram gives r+ — )• ^^^~fi~^ [39] or t ^>- firj [lO] . 
and they are also proportional to tan^ (5/m\. Therefore, the correlation between the Higgs 
mediated r decay and Bg — t- ji^ ji^ can be important to investigate the quark and lepton 
unification. 
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Figure 4: Br(r — )■ yU7) is shown as a function of Br(r — t- /ir/) for universal boundary condition 
(given in the text). The red dotted points satisfy the CDF allowed range of Bt{Bs -^ fi^fi~). 



The current experimental bounds are Br(r — )■ ^rj) < 2.3 x 10~^ [H] and Br(r -^ 3fi) < 
2.1 X 10~® [12]. Because the Higgs mediated contribution to these processes generates a fixed 
ratio Br(r — )■ ;U?7)/Br(r -> 3/u) = 8.4 |40], the r — )■ /i?7 decay gives more important bounds to 
study the Higgs mediated t-^ flavor violation [i3] . 

The measured Bt{Bs — ?■ fJ'^fi~) has an impact on the Br(r — )■ /ir/) for given boundary 
conditions. In FigJH we plot Br(r — )■ firj) and Br(r — )■ yU7) for r7ii/2 = 500 GeV, Aq = 
and tan/3 = 50. We vary the SUSY breaking sfermion mass niQ < 1.5 TeV. We take the 
non- universal SUSY breaking Higgs masses to vary /i and rriA- We assume that the left-handed 
slepton mass matrix can have off-diagonal elements due to the right-handed neutrino loop. Since 
the T-fi Higgs coupling is induced by the finite correction, the t ^ fij operator is correlated 
and they are proportional to each other for given mass spectrum. The red points in the plot 
satisfy the 1 sigma range of Bt{Bs — )■ /^"""/i") reported by the CDF. The CDF result constrains 
Br(r — !■ fit]) for a given value of Br(r — )■ fi'y). 

The Higgs mediated Br-sl coupling is generated by stop-chargino loop diagram, and the 
Higgs mediated tr-^l coupling is generated by stau-chargino loop diagram. As a result, the stop 
and stau mass spectrum is important for the correlation between Br(r — )■ ^rf) and Br(r — )■ yU7) 
for a given Br(i?s — )■ fi~^fi~). Therefore, the correlation can be a probe of the SUSY breaking 
sfermion mass universality and the gaugino mass unification. For example, if the sfermion 
masses are not universal, the correlation between Br(r — )■ firj) and Br(r — )■ fi'j) will be broken. 
In order to illustrate it, we show the plot for the case of sfermion mass non- universality. In 
FiglHl we plot Br(r — )■ fir]) and Br(r -^ ^'-f) for mi/2 = 500 GeV, Aq = and tan/3 = 50. 
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Figure 5: Br(r — )■ fi'-f) vs Br(r — )■ /ir^) is plotted for non-universal boundary condition (given 
in the text). The points satisfy the CDF allowed range of Bt{Bs — )■ fi^fi~). The red dotted 
points satisfy the universal boundary conditions (as in Fig. 4). 



We choose the SUSY breaking sferniion mass rriiQ ^ m^ where m 
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mi. The points satisfies the 1 sigma Bt{Bs — )■ fj^^fi") range reported by CDF. 
The red points correspond to the case of sferniion universality and thus, they are same as the 
red points in FigJH The green points correspond to the case of non-universal sfermion mass. 

The correlation between Br(r — )■ /xr/) and Br(r — )■ fi'j) can be also broken if b-s Higgs 
coupling is accidentally canceled by a new left-handed squark flavor violating source (which is 
absent in the minimal type of SU(5)). If the sizes of the off-diagonal elements are same for all 
sfermion species, the r — )■ /i7 bound can be easily saturated and the Bg — )■ fj^^fi" amplitude will 
not be canceled. Therefore, the flavor violation in the quark sector has to be larger than one 
in the lepton sector choosing the SO(IO) breaking vacua to allow the cancellation. However, if 
the bn-SL Higgs coupling is canceled (namely the new FCNC contribution is destructive), then 
the new FCNC provides constructive SUSY contribution to the 6 — )■ 57 operator (C72,) [15], 
which is disfavored from the experimental constraint. 

It is expected that the upper bounds of the branching ratios of flavor violating r decays can 
be roughly one order below at the super B factory [Hj. It may be hard to achieve the allowed 
region of Br(r — )■ firj). If, however, Br(r — )■ fir]) is measured to be ^ 10~^, the slepton and 
squark mass universality or gaugino mass unification may need to be broken to explain. 



5 Conclusion 

The b-s flavor changing transition is one of the important probes of new physics. The purely 
leptonic Bg meson decays are helicity suppressed in SM, and the prediction of the branching 
ratio is very small. Therefore, the measurement of the rare decay process of Bg — > yU+/i~ plays 
an important role in finding the signature of the new physics beyond SM. In fact, in the MSSM, 
the decay rate has a strong dependence on tan/3, and it can be large even if the colored SUSY 
particles are heavy, as preferred by the current experimental constraints. In other words, a 
large tan /3 is preferred to reproduce the large branching ratio of Bg — )> fi^fi^ for heavy colored 
SUSY particles. In this way, the CDF measurement of the Bg — ?■ fi'^fi~ rare decay process has 
a tremendous impact to bound the SUSY parameters, and can provide a crucial upper bound. 

Even in the minimal models of the SUSY breaking parameters in which flavor universality is 
assumed, the Bg — )■ fi~^fi~ process can be large due to the flavor violation originated from CKM 
mixings. In GUT models, the flavor violation can be related to the symmetry breaking and the 
GUT particle spectrum, and the branching ratio of Bg — )■ fi~^fi~ process can give a constraint 
on the GUT models. Especially for the minimal type of SU(5) GUT where the flavor violation 
comes from the Dirac neutrino Yukawa coupling, the parameter space is really restricted and 
the branching ratios Br{Bg — )■ fJ^^fi") and BR(r — )■ yU7) are bounded from below if the CP 
violation in Bg mixing is large, which is indicated by the Bg — )■ J/ip(f) decay reported by the 
Tevatron and the LHCb, and the dimuon asymmetry in the semileptonic B decays reported by 
DO. The flavor violating lepton decays are also important to flnd a footprint of the GUT scale 
physics. 

The excess of the Bg — )■ fi^fi" decay deviation from the SM prediction at more than 90% 
CL can be soon verifled at the LHC, and the squarks and the gluino can be soon found if the 
true value of Br{Bg — )■ /i'^/U") lies in the 1 sigma range of CDF measurements. The large Bg 
CP phase can be also soon verifled at the LHC. A surge of experimental results will appear 
soon to test the whole structure of the models presented. 
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